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The theory of current transport in a narrow superconducting channel accounting for thermal fluc- 
tuations is developed. These fluctuations result in the appearance of small but finite dissipation 
in the sample. The value of corresponding voltage is found as the function of temperature (close 
to transition temperature T — T c <C T c ) and bias current J < J c ( J c is a value of critical cur- 
rent calculated in the framework of the BCS approximation, neglecting thermal fluctuations). It is 
demonstrated that the value of the activation energy 5F (exponential factor in the Arrenius law) 
when current approaches to the critical one is proportional to (1 — J/J c ) 5 ^ 4 ■ This result is in con- 
cordance with the one for the affine phenomenon of the Josephson current decay due to the thermal 
phase fluctuations, where the activation energy SFj ~ (1 — J/J c ) 3,/2 (the difference in the exponents 
is related to the additional current dependence of the order parameter). Found dependence of the 
activation energy on current explains the enormous discrepancy between the theoretically predicted 
in Ref. 1 and the experimentally observed broadening of the resistive transition. 
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For the first time the role of fluctuations in the en- 
ergy dissipation in the process of current flow through 
the narrow superconducting channel (NSC) was consid- 
ered in the paper of Langer and Ambegaokar— more than 
forty years ago. Publication of this paper has strongly 
influenced all further research in this field, it became clas- 
sical, and corresponding results were included in multiple 
monographs and handbooks on superconductivity^^. 
Nevertheless, even the authors of Ref pi themselves men- 
tioned the striking discrepancy between the predicted 
and the experimentally observed^ values for the width 
of resistive transition. They attributed such discrepancy 
to the possible presence of inhomogeneities in the sam- 
ples. 

Below we will show that point is the Ref. 1 contains 
two incorrect assumptions which result in the paramet- 
rically large overestimation of the activation energy in 
the exponent of Arrenius law. The first one is related to 
the choice of the form of the free energy functional F s , 
where side by side with the standard Ginzburg-Landau 
(GL) part the current-field interaction term should be 
taken into account. The second is related to definition 
of the saddle point in the Arrenius law. The latter 
should correspond to the second stationary solution of 
the GL equation with fixed value of the flowing current 
J, while the authors of Ref i just accepted it in the form 
A(x) = tanh \x/ (v2£gl (T))] , which is correct only 
in the absence of current (£gl (T) is the GL coherence 
length) . 

In this Letter we have calculated the value of activa- 
tion energy 6F for the NSC biased by current J. In 
order to do this we wrote the free energy functional in- 
cluding both GL and the current-field interaction terms, 
derived corresponding GL equations, and found the order 
parameter A(%, J). Thus we will show that taking into 



the flowing current results in the considerable decrease of 
the value of activation energy with respect to the result 
of Refi. For realistic currents this decrease can reach up 
to two orders of magnitude. 

Generalities and stability problem. Let us start our 
discussion considering the free energy functional written 
for NSC biased by current J: 
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Here v = mpp/ (2ir 2 h 3 ) is the density of states (pp is the 
electron Fermi momentum), r = 1 — T/T c is the reduced 
temperature, cL = d/dr—2ieA/c, ( (x) is the Riemann 
zeta- function, c is the speed of light, 
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is the diffusion coefficient (Z tr and r tr are the elec- 
tron transport mean free path and transport scattering 
time)£, ip (x) is the Euler psi-function. In order to avoid 
cumbersome expressions in intermediate calculations we 
will use the system of units where fee = 1 and h = 1. Nev- 
ertheless, in the formulas important for comparison with 
experiment we write these constants explicitly. We as- 
sume the current density j to be constant due to the nar- 
rowness of the channel (its cross-section S <C £q L (T)). 
It is the presence of the bias current J = jS, flowing 
through the superconducting channel, that results in ap- 
pearance of the additional gauge-invariant term in the 
total free energy functional. 

According to general principles, the variation of the 
free energy functional (|TJ) over modulus |A| and gauge- 
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invariant quantity <f> = A 
leads to GL equations 
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Below we will operate in the gauge where A (x) is real 

(v = o). 

The first solution of the system @ with fixed current 
corresponds to the the homogeneous state of the NSC 
with the constant values of the vector potential A and 
the order parameter Aq along the channel: 
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where A 00 (r) = [8tt 2 T 2 t/ (7C(3))] 1/2 is the BCS value 
of superconducting order parameter close to critical tem- 
perature in the absence of current. One can see that the 
current density as the function of vector potential reaches 
its maximal value 
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when the vector potential is equal to |^4 C xtr| 
[2Tc 2 r/ (3ire 2 V)] 1/2 (see Fig. □}. 



For the fixed current value the corrections {Ai, Ai} are 
connected by the simple relation A\ = — 2A(Ai/Ao). 
This relation follows directly from the first equation of 
the system Substituting it to the first equation of 
the system ([2]) one can write the required equation for 
the eigenvalue A : 
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Remaining in the left-hand side value A 2 , can be ex- 
pressed in terms of the vector potential in accordance 
to Eq. (O, which results in 
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The eigenvalue A becomes zero when the vector poten- 
tial reaches its critical value A — A cxtI . This is exactly 
the point of the absolute instability, where the activation 
energy in Arrenius law should turn zero. 

Activation energy in decay rate of the NSC. The system 
© at a given current value has the second, inhomoge- 
neous, solution A(A,x)( x is the coordinate along the 
channel), which determines the value of activation en- 
ergy in Arrenius law. In order to find it let us exclude 
the vector potential from Eqs. ©. One finds 
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FIG. 1: Schematic representation of the current density j in 
the channel as function of the vector potential A. 

Let us study the stability of the found homogeneous 
current state of NSC with respect to growing current. 
As is well known, it is determined by the behavior of the 
eigenvalues of the operator obtained by linearization of 
"equations of motion", in our case, Eqs. ©. The state 
becomes absolutely unstable when the lowest eigenvalue 
turns to zero (see, for example, Rcf. 7 ). In order to find 
corresponding value of the vector potential let us present 
the functions {A, A} in the form {Ao + Ai, A + A\} as- 
suming Ai <C Aq, Ai <C A and linearize the system ©. 
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where C — const. Using the dimensionless variables 

j = j c r, a 2 (r, X ) = a 2 (r) z (x) , a 2 (r) = a 2 £ (r) , 

(6) 



one finds from Eqs. @ the cubic equation for C 
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In the range A < A extI Eq. ([7]) has the only physically 
meaningful solution 
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The value 
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can be found from Eq. ([5]) by applying the boundary 
conditions at infinity: A (oo) = and dA/dxlx^ac = 0. 

Let us note that the solution of Eq. §5§ should be 
even with respect to any fixed point xq : A (x — Xq) = 
A (xq — x) . It is why we can assume Xq = and solve Eq. 
([5]) for x > with a boundary condition dA/dx\ x= o — 
0. Such solution reads as 
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Final integration results in 



Z (x) = — h — tanh 2 | 2x\ 1 
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One can see that corresponding A (r, x) is reduced to 
the one of the Refi only when the flowing current is zero 

(r = o,£ = i). 

Substituting A 2 (T, x) found above to Eq. ([T]) and us- 
ing Eqs. ([H])-© we obtain the expression for activation 
energy SF 
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Simple integration leads to the final expression for acti- 
vation energy, valid for an arbitrary bias current: 



SF = 2vtA 
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When the bias current is close to its critical value one can 
find from Eqs.([9])-([10]) that 3£ - 2 = - T 2 and the 

expression for activation energy is noticeably simplified: 



SF = SF -£ (r) 
where 
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Eq. (jTTj) can be foreseen and seems trivial, while the sec- 
ond cancelation, which results in the additional decrease 
of activation energy SF with respect to SF , is surprising. 

Pre- exponential factor. Let us move to estimation of 
the pre-exponential factor in Arrenius law for the num- 
ber of voltage jumps per unit time. The GL formalism 
does not allow its exact definition: in order to do this 
it is necessary to know at least the dynamical equations 
for the order parameter valid in the wide range of fre- 
quencies. Other possibility is to know the shape of J — V 
characteristic of the NSC above the critical current, for 
J — J c *C J c . Nevertheless, the simplest way to evalu- 
ate the pre-exponential factor is the dimensional analysis 
which we will use below. 

The Josephson relation Ref£ connects the average 
voltage V at the channel to the average time interval At 
between the voltage jumps: eV = ttH/ (At) . The latter 
can be estimated as 
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Indeed, the first factor should define the character- 
istic time scale. We choose it to be in the form 
2nh/kBAf) (r) .Then, one should take into account the 
existence of the "zero-mode" , i.e. the arbitrariness of the 
choice of xq. This means that the instability can arise in 
an arbitrary point of channel and it involves the domain 
of the size of coherence length. This results in appear- 
ance of the second factor in Eq. Q3]), which is nothing 
else as the ratio of the coherence length in the presence 
of current to the length L of the channel. Finally, ac- 
counting for the Arrenius exponent we arrive at the Eq. 

At this point one can write down the J-V characteris- 
tics of the NSC close to transition temperature and for 



is the value of activation energy at zero current (r = 0). 

Let us emphasize that in Refji the activation energy 
in Arrenius law tends to non-zero constant when current 
reaches its critical value. The latter constant differs only 
by the numerical coefficient of the order of one from the 
value of activation energy Eq. (TIB"]) calculated at zero 
current. At the same time, one can clearly see from Eq. 
(Tj"2"]) the dramatic effect on SF of the correct account for 
flowing current. Indeed, the current dependent factor in 
square brackets strongly depletes the activation energy 
Eq. (p~2|) with respect to SFg. Even not too close to the 
critical current when Vl — T 2 « 1/3 the activation en- 
ergy given by Eq. (|12p is 30 times smaller then prediction 
of the Refi (see Fig. |2J. 

Let us indicate the interesting property of the current 
dependent factor in the general Eq. pT|) exposed by 
square brackets. In the vicinity of the critical current 
two first terms of its Taylor expansion are exactly can- 
celed out (see Eq. p^|l ). Cancelation of the first term in 
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FIG. 2: Current dependence of the activation energy. SF (J) 
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arbitrary current J < J c : 
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Discussion.We demonstrated that the account for the 
effect of current flow through the NSC results in a strong 
suppression of the energy barrier for the phase slip events 
with respect to its value at zero current. In the gen- 
eral Eq. not only the first term, proportional to 

1 /4 

(1 — J / J c ) , but also the next one, proportional to 
(1 — J/J c ) 3 ^ 4 are canceled close to the critical current 
J c . As the result, the first non- vanishing term turns out 
to be proportional (1 — J/J c ) 5 ^ 4 , which is the reason of 
a strong reduction of the barrier. Moreover, the addi- 
tional numerical smallness arises due to the high order of 
Taylor expansion in Eq. (jllj) . As a consequence the bar- 
rier reduction turns significant even for currents, being 
relatively far from the critical value: for J c — J = O.IJ c 
the reduction factor is 12.4. 

One can estimate the width of the temperature smear- 
ing AT of the transition at fixed current J just equating 
SF ~ T c . In the most interesting case J — > J c Eq. ([FT 
gives 
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where AT ( J = 0) ~ ^(l) (Ginzburg-Levanyuk number 4 
for the NSC) is the width of transition at zero current. 
This formula differs from the result of Ref. 1 paramet- 
rically, by (1 — J j J c ) 5 ^ 6 , which provides the necessary 
factor of the order of tens lacking in Refi for the agree- 
ment with the experiment 5 -. 

Let us mention that the discussed dissipation process is 
affine to the well studied phenomenon of the Josephson 
current decay due to the thermal phase fluctuations^. 
The activation energy for the latter also turns zero when 
current, flowing through junction approaches the critical 



SFj ~ (1 - J/J c f 2 . 

The difference between exponents of SFj and SF (Eq. 
(fl"2"]) is related to the additional dependence of the order 
parameter in NSC on current. This analogy can be useful 



when one is trying to understand what kind of the J-V 
characteristics one could expect for NSC when currents 
exceed the critical value (J > J c ). Three different scenar- 
ios are possible after overcoming the potential barrier in 
Josephson junction (see Fig. [3]) : (I) the system jumps 
to the neighbor minimum and for a long time remains 
around the new minimum (its phase changes by 2tt); 
(II) the system switches to the regime of the "free over- 
barrier semiclassical motion" ; (III) the system jumps to 
any other minimum with the phase change by 2ttN (N 
is an integer number with some distribution function 10 ) 
and for a long time remains there. 




FIG. 3: Different scenarios for the decay of a metastable state 
in Josephson junction. 



The realization of this or that scenario depends on the 
values of the effective viscosity and J/ J c (depth of the po- 
tential well). One could expect possibility of realization 
of all this variety of options (depending on J/ J C ,T/T C 
and l tr ) also in the supercritical regime of the J-V charac- 
teristics of the NSC. In experimental realization the first 
importance acquires the problem of overheating related 
to the very high current densities in superconductor. 
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